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Вычисление объемных характеристик 
изображений компьютерной томографии” 


В статье предложена иерархическая схема анализа, согласно кот орой вычисление характеристик одного уровня 
проводится на основе результатов предыдущего уровня, а также описан алгоритм вычисления объемных 
характеристик объектов изображений компьютерной томографии на основе связных компонент, ПОЗВОЛИВШИЙ 
уменьшить вычислительные затраты по сравнению с традиционным «объемным» сканированием. 


Введение 


Анализ изображений компьютерной томографии (КТ) [1] так или иначе связан 
с вычислением объемных характеристик, обычно одним из трех способов: на основе 
векторного представления, методом объемной связанной компоненты и путем ап- 
проксимации характеристик плоскостей, образующих объем. Каждый из способов 
имеет свои недостатки. Так, векторное представление З)-объекта всегда имеет погреш- 
ность, вызванную особенностями его реконструкции, метод связной компоненты опери- 
рует растром, что закладывает погрешность, связанную с размером вокселя [2]. Ап- 
проксимация характеристик образующих плоскостей компенсирует этот недостаток, 
предполагает, что между образующими плоскостями объект имеет гладкую поверхность, 
однако медицинские КТ-изображения могут иметь разные межслойные интервалы [3]. 
Вычисления характеристик объекта в бинарном объеме реализуются двумя способами: 
накопления информации во время формирования образа объекта и на основе связанной 
компоненты [4]. 

Целью данной работы является разработка метода вычисления объемных 
характеристик объектов изображений компьютерной томографии в соответствии с 
предложенной иерархической схемой анализа с использованием связных компонент. 

Особенность предложенного алгоритма заключается в организации стека линий так, 
что анализ возможен для изображений с любой связностью. Кроме того, была 
использована иерархическая схема анализа, согласно которой вычисление характеристик 
одного уровня проводится на основе результатов предыдущего уровня. 


Вычисление объемных характеристик первого уровня 
на основе метода связных компонент 


В основе алгоритма лежит идея попиксельного анализа путем растрового ска- 
нирования линий для элементов объекта. Каждая линия при этом сканируется гори- 
зонтально в обоих направлениях от заданного пикселя (рис. 1). 


“ Работа выполнялась в рамках проектов РФФИ 08-01-90016, БРФФИ Т0ЗР-117 и МТНЦ В-1486. 
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Рисунок 1 — Сложный объект и направления его сканирования (прямые линии — 
направления сканирования, штриховые — переход после анализа соседства) 
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Рисунок 2 — Схема алгоритма анализа объекта по связанной компоненте (а); 
схема процедуры анализа линии от заданного пикселя (6) 


Алгоритм. Для входного пикселя сканированием определяется горизонтальная 
линия объекта, и рассчитываются ее координаты, которые заносятся в стек. Одно- 
временно все пиксели линии помечаются флагом или сервисным значением, которое 
исключает их из дальнейшей обработки, и определяются первоначальные характери- 
стики. Затем они берутся из стека и для текущей линии выполняется поиск соседних 
пикселей, принадлежащих объекту. Как только такой пиксель обнаруживается, запуска- 
ется функция анализа и сканирования линии, которой принадлежит найденный пиксель. 
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Новая линия заносится в стек. Для нее рекурсивно выполняется анализ соседства и анализ 
связанных линий. Когда рекурсия заканчивается, из стека извлекаются координаты 
следующей линии и анализ соседних линий повторяется. Алгоритм работает до тех пор, 
пока в стеке не останется линий. В результате получается растровое сканирование 
объекта, программно реализуемое через смещение указателей на пиксели. 

Выше (рис. 2) дана иллюстрация работы алгоритма анализа объекта по связанной 
компоненте и развернутая схема процедуры анализа линии от заданного объекта. 

На схеме (рис. 2 6) КТ - это указатель на крайний пиксель в изображении; 
ВК и ЕВГ - указатели на крайние правый и левый пиксели линии, соответственно; 
у — вертикальная координата линии; хг и х/ — горизонтальная и вертикальная 
координаты правого и левого края линии соответственно, а [епе1й — длина линии. 


Базисные характеристики объектов 


На основе занесенных в стек координат положения линий (у, х,, х,) рассчитыва- 
ются первоначальные базисные характеристики: площадь, периметр, геометрические 
статистические моменты, характеристики яркости и цвета. Эти характеристики 
определяются кумулятивно при сканировании линий. 

Площадь равна сумме всех длин линий на изображении. Периметр — числу 
пикселей на конце линии. Для полноценной линии это 2 пикселя, для одной точки -— 1. 
Нулевой момент и соответствует сумме пикселей внутри объекта и равен традицион- 
ной площади. Моменты первого порядка (и по горизонтали и до: по вертикали), 
нормализированные к площади, соответствуют центру тяжести объекта. Центральными 
называются моменты большего порядка. Они описывают объекты без учета их 
положения в пространстве. На практике используют нормализованные к площади 
центральные моменты, по которым судят о форме объектов. Так, первый центральный 
момент 1 свидетельствует об отличии от круглой формы, вторые центральные 
моменты и/20 И Ио — 0 симметричности объекта относительно горизонтальной и 
вертикальной осей. Для вычисления геометрических моментов используются характерис- 
тики каждой линии объекта: 


ас, = у-1епай; ас, =епай - (х, +х,)/2; ат» =ас, у; ати =4с,-у; 
ат» =1(епей -х, .х, + (а. —епей-х,): (2 {епей —1)/3; (1) 
ат; =ату : уз ату = ас, ас; ат», = ат у; ат» =ат» + у. 


На их основе посредством накопления вычисляются моменты (табл. 2). 
Характеристики яркости определяются как среднее и дисперсия по всем (А, С, В) 
значениям пикселя. 


Таблица 1 — Формулы вычисления некоторых производных характеристик 


Моменты Моменты Моменты третьего 
первого порядка второго порядка порядка 
Но = Мю '@©,; Из = Из ть; Из = Ио ЧТ; 
Ни = Ня -Чс,. ии = И ат; Ив = Из ат; 
Мо = Ид Что. И» = Их `ать; 
И» = Ир `ать.. 
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Производные характеристики объектов 


На втором уровне анализа КТ-изображений алгоритм разветвляется для вычисления 
простых объемных характеристик (объема, поверхности), характеристик, основанных 
на планиметрических моментах (ориентация, анизотропия, оси и т.п.) и для опре- 
деления выпуклой формы фигуры [5], [6]. 

Первая ветвь основана на аппроксимации ряда значений величин (в виде 
неравномерного ряда). Число величин соответствует количеству слоев изображения 
и формируется из значений характеристик плоских бинарных объектов для каждого 
слоя. В случае наличия на одном слое нескольких объектов их значения суммируются. 
К таким характеристикам относятся: объем (на основе площадей слоев), поверхность 
(на основе периметров), выпуклый объем (на основе выпуклых площадей слоев), 
выпуклая поверхность (на основе выпуклых периметров) и т.д. 

Характеристики ориентации объекта и оси [2] вычисляются на основе сумм 
геометрических моментов с учетом ширины между слоями: 


Апё[е = У, апё[е„ Чл › (2) 
п 
где и — номер слоя, Аизе — ориентация объекта, апз{е„ — ориентация объекта на слое, 
4, -— толщина слоя. 

При определении выпуклой формы объекта затруднение представляют слои с 
несколькими объектами, которые связаны друг с другом через другие слои. Решение этой 
проблемы было найдено в объединении всех объектов в один линией, проведенной через 
центры тяжести объектов. В результате строится выпуклый контур на слое, включающем 
все объекты. Площадь поверхности, площадь выпуклой поверхности, объем и выпуклый 
объем находятся интегрированием периметров, выпуклых периметров, площадей и 
выпуклых площадей соответственно. На основе полученных объемных характеристик на 
третьем уровне вычисляются значения таких комплексных параметров (табл. 2), как: 

— компактность (сотрасте$$, или сте) как отношение объема объекта к объему 
шара с одинаковой площадью поверхности. Шар считается объектом самой компактной 
формы, для него компактность равна 1; 

— сферичность (5рйетсйу), свидетельствующая о локальной нерегулярности. Признак 
сферичности вычисляется как отношение объема объекта к объему описанной сферы, 
равной по площади поверхности объекта; 

— выпуклость (сопуехйу) как отношение выпуклой поверхности к обычной. Несмотря на 
большое количество других известных характеристик, признак выпуклости является 
достаточно информативным при описании формы объекта. Если объект выпуклый, эта 
величина будет равна 1, при появлении впадин она уменьшается; 

— целостность (5оПаЙу) — вторичный признак, описывающий структурную плотность и в 
некотором роде однородность объекта. Он определяется как отношение обычного и 
выпуклого объема. Если объект целый и не имеет дыр и впадин, эта величина равна 1, 
при их наличии она уменьшается. 


Таблица 2 — Формулы вычисления некоторых производных характеристик 


компактность сферичность выпуклость целостность диметр 
экв. шара 
2.л’’?.ГОГИМЕ 21°? .ГОГОМЕ сопуех зи {асе ГОГИМЕ 3 ГОРИМЕ 
3. зитГасе'' 2 3. (сопуехзи асе!" зи асе сопуех уоште 2.л 
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Дополнительно определяются объемные координаты объекта как минимальные и 
максимальные значения. Усреднением значений для плоскостей определяются центр 
масс, средняя яркость и оптическая плотность объекта. Существует еще ряд полез- 
ных производных характеристик, например, диаметр эквивалентного объекту шара. 


Выводы 


Иерархическая схема анализа КТ-изображений дала комплексное решение 
задачи вычисления более 90 характеристик объектов разной сложности по восьми и 
четырехсвязности. Модификация, внесенная в алгоритм послойного определения 
объекта по связанной компоненте через стек, позволила увеличить скорость работы 
в 1,76 раза по сравнению с простым объемным сканированием. 
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